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ABSTRACT 

Exoplanets which are detected using the radial velocity technique have a well- 
known ambiguity of their true mass, caused by the unknown inclination of the plan- 
etary orbit with respect to the plane of the sky. Constraints on the inclination are 
aided by astrometric follow-up in rare cases or, in ideal situations, through subse- 
quent detection of a planetary transit. As the predicted inclination decreases, the 
mass of the companion increases leading to a change in the predicted properties. Here 
we investigate the changes in the mass, radius, and atmospheric properties as the in- 
clination pushes the companion from the planetary into the brown dwarf and finally 
low-mass star regimes. We determine the resulting detectable photometric signatures 
in the predicted phase variation as the companion changes properties and becomes 
self-luminous. We apply this to the HD 114762 and HD 162020 systems for which 
the minimum masses of the known companions places them at the deuterium-burning 
limit. 
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1 INTRODUCTION 

' Planets discovered using the radial velocity technique con- 
, tinue to form a major component of the known exoplan- 
1 ets. These planets have a well-known ambiguity to their 
' masses due to the unknown inclination of their orbits. 
, Some planets have subsequently been found to be more 
massive than originally thought when constraints are later 
placed upon their inclinations. These constraints can come 
from dy namical considera tions, such as for the HD 10180 
system llLovis et al.l 20l"ll), t hrough astrometric follow-up 
l|RefFert fc Quirrenbachl 201 1 | ). or measurements of the pro- 
jected equatorial velocity I Watson et al.l |2010| ). although 
some hot Jupiters have been found to exh ibit spin-orbit mis- 
alignment (see for example IWinn et all (j200^)). The con- 
squence of these constraints can be either to confirm their 
planetary candidacy or to move the mass into the regime of 
brown dwarfs and low-mass stars. 

The mass-radius relationship of short-period exoplan- 
ets is evolving t hrough the discovery of numerous tran- 
siting exoplanets (Burrows et al. 1 12007^ : iFortnev et allbOOTl : 
ISeager et al.l 12007 V The understanding of this relation- 
ship is undergoing continued and rapid evolution through 
the r elease of transiting planets from the Kepler mis- 
sion (|Borucki et al.l l2011al lbl). Surveys for transiting exo- 
planets have also led to the serendipitous discovery of 
transiting brown dwarfs, such as CoRoT-3b (jDeleuil et al.l 



l2008l), WASP-30b ^Anderson et all [201 ll ). and LHS6343C 



( Johnson et ablbOllT ). The mass-radius relationship of low- 
mass stars ha s been investigated by nume r ous authors 
dDemorv et a l.l l2009l : iFernandez et al.l l2009l : iKraus et al.l 
120111 : IRibasll200d ) through new radii measurements and the 
development of models to explain the observed relationship. 
Even so, the low-mass stars for which accurate radii have 
been determined remains a relatively small sample from 
which to draw upon to derive the theoretical framework for 
model construction. These mass-radii distributions give us 
clues about the formation mechanisms that occurred within 
these systems. 

The phase functions and resulting photometric light 
curves of orbiting companions depends upon their 
radii and albedo as well as the orb it al components 
(iKane fc Gelino' '2010"; 'Seag er et aP l2000l : ISudarskv et^ 
12005^ . Kane & Gelino (201 laj') showed how the phase curves 
of exoplanets varies with inclination. However, it was as- 
sumed that the fundamental physical properties of the 
planet (such as mass, radius, and atmospheric properties) 
remain the same. This is particularly relevant to planets 
discovered using the radial velocity technique since these 
planets are inherently subjected to an ambiguity in the mass 
with respect to the unknown orbital inclination. Here we ex- 
pand upon this topic to investigate how the properties of a 
known companion can be expected to change as the inch- 
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nation is decreased and therefore the mass is increased. If 
the companions are self-luminous then this places an extra 
constraint upon what one can expect to see in high-precision 
observations designed to measure phase variations. We thus 
produce a criteria from which one can distinguish between 
stellar and planetary companions solely from high-precision 
photometric monitoring and without the need for astromet- 
ric observations. 



2 INCREASING THE COMPANION MASS 

In this section, we describe the changing properties of a 
companion as the mass increases from the planetary regime, 
through the brown dwarf regime, and into the realm of low- 
mass stars. We further investigate how these properties in- 
fluence the photometric properties up to the self-luminous 
threshold. 

2.1 Mass-Radius Dependence 

The mass-radius relationship of exoplanets, brown dwarfs, 
and low-mass stars is a difficult subject to relate to mea- 
surements since these require a system where the orbital in- 
clination allows for the observations of eclipses and transits. 
Locating such systems involves the monitoring of many stars 
since the probability of having a favourable orbital inclina- 
tion tends to be relatively low, depending upon the semi- 
major axis of the orbiting companion. 

The radii for objects of a given mass depends upon a 
number of factors, including age and metallicity. Here we 
consider host stars which comprise the bulk of the known 
radial vel ocity hosts which h ave typical ages in the range 
1-5 Gyrs (Takcda ot al. 2007). A study of 49 exoplanet host 
stars by Saffc ct al. (2005) found median ages of 5.2 and 
7.4 Gyrs, using chromospheric and isochrone methods re- 
spectively, with dispersions of ~ 4 Gyrs. We thus consider 
objects which are of order 5 Gyrs old. 

The true mass of companions detected using the ra- 
dial velocity technique is the measured mass divided by sin i 
where i is the orbital inclination relative to the plane of 
the sky. As the inclination decreases from an edge-on con- 
figuration (i — 90°), the true mass increases which in turn 
effects the radius and atmospheric properties of the compan- 
ion. It is generally held that the separation between planets 
and brown dwarfs is the onset of deuterium burning. How- 
ever, the mass criteria for determinaing where deuterium 
burning occurs can be quite broad depend ing upon he- 
lium I deuterium abundances and metallicity (|Spiegel et al.l 

By way of demonstration. Figure [T] plots the C OND 
evolutionary model isochrones of iBaraffe et al. I l|2003l ) for 1, 
5, and 10 Gyrs and for a mass range of 0.003-0.100 Mq. Also 
shown are the masses of two known exopl anets, HP 114762 b 
l|Latham et al.lll989l ) and HD 162020b l|Udrv et al.ll2002l ). 
Assuming an inclination of i = 90°, the measured mass of 
HD 114762b is 11.0 Mj (0.0105 Mq) whereas for a near 
face-on inclination of i = 10° the mass increases to 63.2 Mj 
(0.0604 Mq). According to the 5 Gyr isochrone, this re- 
duces the predic ted radius from 0.101 Rq to 0.079 Rq. 
iKane et al.ll|201ll ) have excluded transits for this companion 
which thereby restricts the inclination to be less than 89°. 
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Figure 1. The mass-radius relationshi p of planets t hroug h to 
low-mass stars using the isochrones of iBaraffe et al.l ||2003|) for 
1, 5, and 10 Gyrs. The increasing masses of HD 114762b and 
HD 162020b are shown to demonstrate the effect of decreasing 
the inchnation from an edge-on orbit (i = 90°) to an almost face- 
on orbit {i = 10°). 

For HD 162020b, the same change in inclination increases 
the measured mass of 14.4 Mj (0.0137 Mq) to 82.9 Mj 
(0.0792 Mq), thus decreasing the predicted radius from 
0.099 Rq to 0.097 Rq. In each case, the change in incli- 
nation moves them from the planetary regime to the brown 
dwarf/low-mass star boundary and beyond. 

As noted in Section[T] our knowledge of companion radii 
in different mass regimes is a currently evolving topic. For 
example, an important consid eration is the existence of the 
so-called brown d warf desert (jGrether fc Lineweaveij I2OO6I : 
iKraus et aLl B) which may bias inclinations away from 
this mass regime. The COND evolutionary models utilized 
above will undoubtedly undergo slight adjustments as the 
sample size of known low-mass transiting companions in- 
creases. One of the important parameters is the age of the 
host star. As shown in Figure [1] the range of radii can be 
quite diverse for a given mass depending on which isochrone 
one adopts for the host star. However, the divergence is most 
significant for relatively young (~ 1 Gyr) stars beyond which 
the radii converge upon a small range of radii for a large 
range of ages. Even so, one must consider the uncertainty 
in the host star age, assuming that the companion is of a 
similar age. 

2.2 Impact on Phase Curves 

For companions which are not self-luminous, there will be 
a photometric phase signature from the companion whose 
amplitude depends upon a variety of factors including the 
companions radius, semi-major axis, eccentricity, orbital in- 
clination, and atmospheric pro perties (geometr i c alb edo). 
Here we adopt t he form alism of iKane fc Gelind l|2010l ) and 
iKane fc Gelind (l2011al ) to demonstrate the impact of the 
changing companion properties on phase curves. 
The phase angle of the planet (a) is given by 

cos Q = sin(cj -I- /) sin i (1) 

where lo is the argument of periastron and / is the true 
anomaly. This angle is defined to be a = 0° when the planet 
is on the opposite side of the star from the observer. Note 
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that this means that cos a = 1 is only possible for edge-on 
orbits. All other inclinations will result in a more compli- 
cated dependence on th e orbital parameters as described by 
iKane fc Geiiiiol lloili). 

The flux of a planet fp and host star /* has a ratio 
defined as 



(2) 



/*(A) 

where the flux is measured at wavelength A, Ag{X) is the 
geometric albedo, g{a, A) is the phase function, and Rp is 
the radius of the planet. The star-planet separation r is 
given by 



a(l-e^) 
1 + e cos / 



(3) 



where a is the semi-major axis and e is the orbital ec- 
centricity. The phase function is primarily dependent upon 
the phase angle, whether one assumes a Lambert sphere 
or something more complicated to describe the scatter- 
ing properties. The geometric albedo for gas giant plan- 
ets depends upon the incident flux (star-planet separa- 
tion) since this can determine the amount of reflective 
condensates that can maintain a presence in the up per 
atmospheres ljKane fc Gelincill201ol : ISudarskv" et al.ll2005l ). 

A component of Equation[2]which explicitly relies upon 
the companion properties is the radius squared. As shown 
in Section [2. 11 variation of the companion mass in the range 
0.003-0.100 Mq can vary the radius by as much as ~ 20%. 
This results in a variation of the flux ratio amplitude by as 
much as ~ 40%. The uncertainty in both the planetary ra- 
dius and albedo are discussed further in Section [G] For now 
we note that there is an apparent degeneracy in the pre- 
dicted radius as a function of mass in the isochrones shown 
in Figure [T] due to the minimum located at 0.06 Rq. How- 
ever, beyond this value the companions start to become self- 
luminous. This is discussed in further detail in Section [3l 



2.3 Ellipsoidal Variations 

An additional possible effect of an increased mass is the 
induct ion of ellipsoidal variations in the host star. iDrakd 
l|2003l ) have proposed the use of ellipsoidal variations as a 
tool in eliminating false positives due to eclipsing binaries 
from transit surveys. Stellar ellipsoidal variations have been 
previously detected for short-period giant planet systems 
though, such as the case of HAT-P-7b which was observed 
to have sign atures of pha s e as w ell as ellipsoidal variations in 
Kepler data lWelsh et al.l (|2010l ). The detection of ellipsoidal 
variations i n the contex t of th e Kepler mission is discussed 
in detafl bv lPfahl etal] (|2008h . 

Comparison of phase and ellipsoidal variation am- 
plitudes to those induced by reflex Doppler motion has 
been undertaken in sev eral stud ies, including those by 
iLoeb fc GaudU (|2003t ) and lZucker et all (|2007l '). A full expan- 
sion of the discrete Fourier series for the ellipso i dal va riation 
amplitude is provided by iMorris fc NaftilarJ J 19931). Here 
we ad opt the approximate relation used bv iLoeb fc Gaudil 
lj2003t ) as follows 
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Figure 2. The mass dependence of the predicted ellipsoidal vari- 
ation amplitude and the same range of masses as shown in Figure 
[T]for three different periods. 



where /3 is the gravity darkening exponent. For the purposes 
of demonstrati on, we have assumed /3 = 0.32 as determined 
bv lLucvl (Il967l '). We refer the reader to IClaretl (2000 ) for a 
much more thorough treatment of the dependence of /3 on 
various stellar properties. Figure [2] shows how the ellipsoidal 
variation amplitude can be expected to change as a function 
of companion mass for three different orbital periods. The 
scale used on the x-axis is identical to that in Figure [T] As 
one would expect from Equation |4l a change in the compan- 
ion mass results in an almost linear change in the observed 
flux ratio changes ("almost" since changing the mass also 
changes the semi-major axis at which the expected period 
occurs). Changing the period has a much more dramatic ef- 
fect. Note however that a low-mass star in a 10 day orbit 
will induce a similar ellipsoidal variation to a giant planet 
in a 5 day orbit. 

There are numerous examples of a changing tidal dis- 
tortion a mplitude due to eccentr ic orbits, such as those con- 
tained in lSoszvnski et al] l|2004 ). To account for orbital ec- 
centricity and inclination, we modify Equation |4] as follows: 



^mAt) H'(^ + /)+sin^(^ + /)cos2i)^ (5) 



Fo ^M* V a J 



(4) 



AF 

where we have replaced a with r, as per Equation^ If there 
is signiflcant orbital eccentricity present, one can use this 
as an additional diagostic in interpreting the uniqueness of 
the ellipsoidal variation solution. However, eccentric orbits 
are far more likely to occur at longer periods where the 
amplitude of the variation is greatly reduced. 



3 CROSSING THE SELF-LUMINOUS 
THRESHOLD 

Moving from the brown dwarf into the low-mass star regime 
introduces many changes associated with the companion be- 
coming self-luminous. Planets and brown dwarfs are typi- 
cally not self-luminous except at a young age (< 1 Gyr). 
For brown dwarfs engaged in deuterium burning, the transi- 
tion to hydrogen burning will abort if the mass is less than 
~0.07 Mq after which the object will undergo luminosity 
decay beyond an age of ~1 Gyr. A full description of object 
characteristics in terms of age, helium/deuterium fraction. 
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and metallicity can be found in iBurrows et ai] (|l993l ). Here 
we briefly discuss the main features that one can expect to 
observe as the object crosses the self-luminous threshold. 



3.1 Elimination of the Phase Function 

A major feature of planetary phase functions is that it relies 
on reflected light which is incident from the host star. If 
the dominant source of light from the companion becomes 
that which originates from self-luminosity then there will 
no longer be a phase function of the kind described earlier. 
To account for this, we include a gradual decrease in the 
phase variation amplitude for companion masses > 50 Mj. 
There will however be an increasing likelihood of observable 
ellipsoidal variations with increasing mass (see Section [2. 3p 
such that these variations may subsume the decreasing phase 
effects. 

An additional effect which may produce a phase- like sig- 
nature is the "reflection effect"; a well-known phenomenon 
for close binary stars in which the irradiation of one compo- 
nent by the other r esults in a diff erential flux over the irradi- 
ated stars surface (jWilsonll 19901 ). Application of this process 
to binary systems consisting of a star and a g iant planet 
or brown dwarf is considered by iBudail (|201lh . although 
iGreen et al.1 (|2003l ) argue that this heating of the primary 
star due to the presence of even a close-in gas giant planet 
will be negligibly small. The irradiation will only become 
comparable to the amplitude of the phase variation when 
the mass enters the low-mass star regime, at which point 
the convective atmosphere will in turn ha ve an increased 
albed o due to irradiation from the primary (|Harrison et al.l 
I2OO3I I. 

For planetary and brown dwarf companions, the heat- 
ing of the day-side may result in signiflcant thermal flux be- 
ing re-radiated. This depends not only on the incident flux 
but also the fundamental atmospheric p roperties which gov- 
ern radiative and advective time scales |Fortnev ct al. 2003; 
iKane fc Gelindl2011bl ; IShowman et al.ll2009l )." Thus cfflcient 
atmospheric re-circulation of the received heat to the night- 
side of the planet may be sufficient to suppress an overall 
phase function at infra-red wavelengths. 



3.2 Introduction of Stellar Activity 

At the level of photometric precision required here, it is im- 
portant to consider the level of intrinsic stellar variability. 
Depending upon the properties of the companion, there may 
be intrinsic stellar activity which reveals its nature through 
precision photometry. There are now hundreds of known 
brown dwarfs with an in creasing understanding of their at- 
mospheres and activity ijKirkpatrick et aDl201ll ). However, 
activity and rotation rates amongst brown dwarfs and low- 
mass stars can vary greatly, mostly due to the large range 
of temperatures and ages of those objects monitored. Sur- 
veys of relatively young low-mass stars tend to show signs 
of high chromospheric magnetic activity an d increa sed rota- 
tion rates with dec r easing temperatur e ([Becker et al . 2011; 
Ijenkins et al.|[2009l : IScholz et al. 2011). A spec tral analysis 
of an M dwarf sample bv lBrowning et al] l|2010l ) showed that 
increased rotation appears to be more common in stars later 
than M3, indicating that measurable rotational braking is 



reduced for fully convective stars. iBerta et al.l (|201lh inves- 
tigated the variability of the exoplanet hosting star GJ 1214 
(0.16 M0) and find a periodicity of 53 days with an am- 
plitude of 3.5 millimags at the MEarth bandpass. Clearly if 
one has the precision to detect planetary phase variations 
then one also has the capability to detect stellar activity 
from the companion. Indeed it is expected that a low-mass 
star instead of planetary companion will show larger chro- 
mospheric a nd coronal activity com pared with equivalent 
single stars jZagarashvili et aLlBoO^ ). If the orbital param- 
eters of the companion are sufficiently understood from ra- 
dial velocity measurements then one can extract the vari- 
ability signals of the companion and the host star. There 
may be cases where the orbital period of the companion 
is close to the rotation period of the star, which is gener- 
ally in the range 10- 40 days for radial velocity host stars 
(Sim pson et al.l2010l '). In such cases, the peaks in the power 



spectrum from a fourier analysis of the photometry may sep- 
arate to a degree where starspot variability due to rotation 
can be isolated. 

This discussion is referring to the activity of a potential 
stellar-mass companion. However, the intrinsic stellar vari- 
ability could well substantia lly stronger than th ese effects. 
An analysis of Kepler data bv lCiardi et al.l l|201ll) found that 
most dwarf stars are stable down to the the precision of the 
Kepler spacecraft, with G-dwarfs being the most stable of 
the studied spectral types. The main cause of photometric 
variability in F-G-K stars is starspots and rotat ion, as veri- 
fied by the Kepler variability study performed bv lBasri et al.l 
(|2011 ). Since the orbital period of the radial velocity target 
stars is well determined, this will aid in separating the sig- 
nals of planetary phase from that of the host star variability 
whose period is likely related to the stellar rotation period. 
It should be noted that disentangling the variability of the 
host star may result in an increased the observing time re- 
quirement. 



3.3 Doppler Boosting 

Doppler boosting occurs due to a relativistic effect which 
creates an apparent increase and decrease in the light from 
the host star when it is moving towards and away from the 
observer. The fractional amplitude of the effect is given by 



(6) 



Af _ (3 - a)K 
'Fh ^ c 

where K is the radial velocity semi-amplitude and a is 
the derivative of the bolometric flux with respect to the 
frequency in the stationary frame of reference. This ef- 
fect has been discussed in terms of planetary compan ions 
and high-precision photometry bv lLoeb fc Gaudl (|2003l ) and 
Faigler & Ma zeh (2011,). Doppler boosting due to stellar- 
mass companions ha s been detected in the photometry from 
the K epler mission l|Faigler et al.ll2012l :l van Kerkwiik et al.l 
|20l3). The photometric variability, including Doppler boost- 
ing, induced by KO I-13.01 led to it s conflrmation as a planet 
(Mazch ct al. 2012: | Shporer et a"l] |201 1). It was addit i onally 
shown by iLoeb fc Gaudil (|2003l ) and Ishporer et al.l (l201lh 
that the effect of Doppler boosting can be significantly larger 
than those from phase and ellipsoidal variations for orbital 
periods greater than ~ 10 days. 
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Figure 3. Light curves from Quarter 2 Kepler data for two cases of ellipsoidal variations, Kepler IDs 8386198 (left panel) and 9544350 
(right panel). The fluxes have been converted to perecentage variations against the mean flux of the data. 



We are considering planets which have been detected 
with precision radial velocities. The amplitude of Doppler 
boosting is directly proportional to K and thus we remove 
this effect with only one free parameter. Note that observ- 
ing this effect does not remove the sini ambiguity of the 
companion mass. This is discussed further in Section [l] 



3.4 The O'Connell Effect 

The O'Connell effect refers to the height difference that may 
occur between the maxima in the hght curves of close binary 
stars. There are several possible causes for this effect, such 
as starspots or gas streaming between the binary compo- 
nents. This eff ect has been previously s tudied for eclipsing 
bi narv stars bv|Davidge fc Milond (|l984l ) and more recently 
by IWilsev fc Beakvl l|2009i r Binaries with low-mass stellar 
comp anions are known to exhibit this effect (|Austin et al.l 
I2OO7I ). making this a noteworthy phenomenon here. In par- 
ticular, the precision of data from the Kepler mission allows 
an unprecedented investigation of the frequency and source 
of this effect. 

In Figure [3] we show two example hght curves which 
show several of the effects discussed in this section, Kepler 
IDs 8386198 and 9544350. These data are corrected flux val- 
ues from the Quarter 2 release of the Kepler mission, ex- 
tracted using the Kepler interface of the NASA Exoplanet 
ArchiveQ- These stars were included in t he Kepler Eclips- 
ing Binary Catalog by IPrsa et al] (|201ll ) as having signa- 
tures due to ellipsoidal variations but they can also be seen 
to have clear O'Connell effect signatures. The amplitude of 
these variations are substantial: 5-10 millimags, and so are 
well within detection thresholds for stellar-mass companions 
(see Figure [2]). This is also significantly larger than the pre- 
viously mentioned variations and so care must be taken to 
distinguish this effect for those more subtle in amplitude. 



^ http:/ /exoplanetarchive. ipac.caltech.edu/ 

^ http://astro4.ast.villanova.edu/aprsa/kcpler/ 



4 OBSERVABLE SIGNATURES 

Here we combine the effects we have considered in previ- 
ous sections and apply them to the two examples discussed 
earlier, HD 114762b and HD 162020b. Recall that these var- 
ious observable signatures are being considered for compan- 
ions which are known through their precision radial veloc- 
ity detections. In the case of the two examples being con- 
sidered here, the radial velocity amplitudes are 612 ms~^ 
and 1813 ms~^ which leads to fractional Doppler boosting 
amplitudes of 8.8 x lO"'' and 2.6 x 10~^ respectively. For 
the orbits of these planets, we assume that Doppler boost- 
ing effects may be well characterized and we inject the two 
remaining and comparable effects of ellipsoidal and phase 
variations. 

It was shown bv lKane fc Gelinol (|2011al ) that the phase 
signature one expects from a planet in an eccentric orbit 
is intricately related to both the inclination and periastron 
argument of the orbit. If one also includes the companion 
mass, radius, and various other photometric signatures as 
free parameters then it can become difficult to disentangle 
the physical model to account for the observations. Shown 
in Figure U is an example sequence of photometric signa- 
tures from a simulated companion in a 5 day orbit around 
a solar-type star. The companion has a minimum mass of 
Mp sin i = 10 Mj and is in an eccentric orbit (e — 0.2) where 
periastron passage occurs between the observer and the star 
{uj = 90°). In each panel, the dashed line is the phase varia- 
tion, the dotted line is the ellipsoidal variation, and the solid 
line is the total combined signature. The flux ratio refers to 
the apparent change in flux of the host star, either by re- 
flected light by the companion or the ellipsoidal distortion 
of the stellar shape. An orbital phase of zero corresponds to 
when the companion is at superior conjunction. Thus the 
phase variation tends to peak near phase zero whereas (for 
edge-on orbits, depicted in the top panel) the ellipsoidal vari- 
ation has two minima: when the companion is either behind 
or in front of the star. At this inclination the companion 
is considered to be of planetary mass and the phase sig- 
nature dominated the photometric variability. The middle 
panel shows the effect of reducing the inclination to 30° such 
that the mass of the companion has now doubled. The phase 
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Figure 4. Photometric signature for a simulated companion in a 5 day eccentric (e = 0.2, u) = 90°) orbit around a solar-type star. 
An orbital phase of zero corresponds to when the companion is at superior conjunction. The dashed line is the phase variation, the 
dotted line is the ellipsoidal variation, and the solid line is the total combined signature. The three panels show the effect of changing 
the inclination of the system from edge-on (i = 90°, top panel) to the smaller inclinations shown in the bottom two panels. 



amplitude is now much reduced and the ellipsoidal distor- 
tion of the star as the companion passes through periastron 
passage is more visible. At an inclination of 10° (shown in 
the bottom panel) the orbit is now observed almost face-on. 
The companion is a high-mass brown dwarf and the phase 
signature is almost eliminated, partly due to the slightly re- 
duced size but mostly due to the reduced contrast between 
the day and night-side of the companion as it becomes more 



self-luminous. Note also the greatly increased amplitude of 
the ellipsoidal variations. 

4.1 HD 114762b 

Shown in Figure [5] are the predicted photometric variations 
expected for the companion orbiting HD 1147 62. We use 
the st ellar and orbital parameters measured bv lKane et al.l 
(|201ll ) . This companion has a minimum mass of Mp sin i = 
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HD 114762 b i = 30.0° M„ = 22.0 M, R„ = 0.9 R, 

p J p J 




Orbital Phase 

Figure 5. Predicted photometric signatures for the companion to HD 114762. The dashed line is the phase variation, the dotted Hne is 
the ellipsoidal variation, and the solid line is the total combined signature. The three panels show? the effect of changing the inclination 
of the system from edge-on (i = 90°, top panel) to the smaller inclinations shown in the bottom two panels. 



10.98 ± 0.09 Mj and is in a ~ 84 day period orbit with an 
eccentricity of 0.34. The periastron argument ofuj — 201.28° 
means that periastron passage takes place almost behind the 
star where the orbital phase is ~ 0.95. Recall from Equations 
[5]and[S]that phase and ellipsoidal variations are proportional 
to and respectively. Thus in this case the phase 
variations dominate the signature until the mass is pushed 
high into the brown dwarf regime as shown in the bottom 
panel. 

There are several aspects of note here. Due to the 



relavitevly large star-planet separation there is only a small 
change in the total flux. However, as the phase signature 
disappears and the ellipsoidal signature grows the total flux 
stays approximately the same. Notice that there is a phase 
offset between the peak of ellipsoidal and phase variations. 
The peak of the ellipsoidal variation depends on the star- 
planet separation and the orientation of the major axis of 
the distorted stellar profile, whereas the phase variation also 
depends upon the phase function of the planet. This is the 
key which unlocks the difference between the planetary and 
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HD 162020 b i = 90.0° M„ = 14.4 M, R = 1.0 R, 

p J p J 




Orbital Phase 

HD 162020 b i = 30.0° M„ = 28.8 M, R„ = 0.9 R, 

p J p J 




Orbital Phase 



Figure 6. Predicted photometric signatures for the companion to HD 162020. The dashed line is the phase variation, the dotted Hne is 
the ellipsoidal variation, and the solid line is the total combined signature. The three panels show the effect of changing the inclination 
of the system from edge-on (i = 90°, top panel) to the smaller inclinations shown in the bottom two panels. 



stellar companion signatures in this case since one can de- 
tect the change in phase far more easily than one can detect 
the change in amplitude of the signature. 

4.2 HD 162020b 

Figure [6] shows the equivalent predicted photometric varia- 
tions for the companion to HD 162020. The minimum mass 
is Mp sin i = 14 A Mj and the orbital period and eccentric- 
ity are 8.43 days and 0.277 respectively with a periastron 



argument oi uj — 28. 4°. These orbital p arameters are from 
the measurements of lUdrv et al.l (|2002i ). who estimate the 
host star as being K3V. We derive a radi us for the host star 
by app lying the surface gravity logy from lValenti fc Fischer! 
l|2005l ) to the relation 

^°Sf' = l°g(^)-21og(^^)+log50 (7) 

where logp© = 4.4374 l|Smallevll2005l ). From this we calcu- 
late a stellar radius of = 0.52 _R0. We additionally cal- 
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culate the radius using the relations of [Torres et al.l (|201Cll ) 
and find = 0.53 ± 0.03 Rq. The eUipsoidal variations 
are sensitive to the radius of the host star (see Equation O 
but the consistency of the radius determinations and small 
uncertainties means that the calculated amplitude is quite 
robust. 

The panels of the figure show that the phase varia- 
tions are able to dominate the total signature even when 
the inclination is decreased to 30°. Ifowever, as the mass 
rapidly increases beyond this point the phase variations drop 
to zero. By the time an inclination of 10° is reached it is 
clear that the companion is able to sustain hydrogen burn- 
ing and the photometric variations and now due to tidal 
distortions induced by the eccentric orbit. As was the case 
with HD 114762b, the total amplitude of the variations does 
not significantly change but the phase offset grows as the el- 
lipsoidal variations emerges as the dominant effect. 



5 A NOTE ON SPECTRAL LINE DETECTION 

An additional effect of increasing the companion mass is the 
introduction of absorption bands which are characteristic 
of brown dwarfs and late-type stars. One may then wonder 
if such features would become apparent in high-precision 
spectral data or if the companion and its star could be con- 
sidered a spectroscopic binary. The optical spectrum of an 
M dwarf is typically dominated by TiO and VO absorp- 
tion bands. Considerable progress has been made on char- 
acterizing the properties of the cooler L dwarfs. An exam- 
ple of an L dwarf was detected using both radial velocities 
and hig h-resolution ima ging orbiting the solar analog IfR 
7672 bv lLiu et all (|2002f ). The L dwarf prototype, GD 165B, 
has been used to model the spectra of many similar types 
of brown dwarfs through tlic identification of prominant 
spectral features (jKirkpatrick et al.lll999i ). In particular, L 
dwar fs exhibit strong metal h ydride bands and alkali metal 
lines. iKirkpatrick et al.l j 19991 ) found the Na doublet at 8183, 
8195 A to be especially strong and a possible candidate 
for detection in optical passbands. For the even cooler T 
dwarfs, such as Gl 229B, the methane absorption bands 
and broad absorption features due to alkali metals tend to 
dominate near-infrared (N IR) spectra l|Saumon et al.ll2000l : 
ISengupta fc Krishanll2000l '). The major hinderance facing de- 
tection prospects for these features is that they primarily 
manifest at red wavelengths where the fiux is dramatically 
reduced at optical passbands. Thus, resolution of the com- 
panion properties through spectral line detection is better 
suited to the near-IR instruments which are being developed 
to perform searches for exoplanets around late-type stars. 



or Doppler Boosting variations. The simulations shown in 
Section|4]demonstrate that the contribution of the planetary 
phase to the total variations become significantly less when 
the orbit is close to face-on, even for eccentric orbits. Thus 
the uncertainties in the radius and albedo become particu- 
larly important for one to be able to discern the planetary 
properties for low inclination orbits. 

Regarding planetary albedo, there has been consider- 
able effort to produce analytic models for estimating the 
albedo at a g i ven st a r-planet separation, such as the w ork of 
ICahov et all (I2OI0I '). iMadhusudhan fc Burro^ (|2012l '). and 
S elsis et al.l l|201lf ). The situation is more complicated than 
a simple distance-albedo relation however, as shown by the 
recent discovery of a suprisingly hi g h alb edo for Kepler- 7b 
iDemorv et al.ll201ll ). ICowan et al.l (|2012l ') provides a sum- 
mary of the recent geometric albedo measurements, the di- 
versity of which may partially be a function of the heat 
redistributi on efficiency in the high equilibrium tempera- 
ture regime. [Pemorv et al.l (|201ll ) attribute the anomalously 
high albedo of Kepler- 7b as being due to a combination of 
Rayleigh scattering and clouds. It is also possible that a 
smooth dependancy of albedo on star-planet separation is 
broken by phase transitions whereby removed reflective con- 
densates reappear in the upper atmosphere as clouds for a 
small range of equilibrium temperatures, thus increasing the 
albedo. Clearly we require a greater understanding of giant 
planet albedos in order to be able to unambiguiously extract 
the phase variations component for non-transiting planets. 

Regarding planetary radius, our understanding of the 
mass-radius relationship is developing rapidly as pointed out 
in Secti on[Tl Based upon the c urrent knowledge of transiting 
planets. iKane fc Gelinol (|2012i ) show that the radii of planets 
with masses greater than ~ 0.3 Jupiter masses follow an ap- 
proximately linear model. This mass range also encompasses 
the majority of non-transiting planets for which their plan- 
etary status is ambiguious. It is expected that the variation 
seen within this region is due in no small part to the fact that 
highly-irradiated giant planets dominate this sample due to 
the bias of the transit technique. However, Kepler is improv- 
ing our understanding of planetary radii at a larger range 
of period than that which is encompassed by the ground- 
based surveys. The proposed Transiting Exoplanet Survey 
Satellite (TESS) will add to this understanding but will also 
beneflt from the follow-up potential of the planet discoveries 
due to the brightness of the host stars (Rickcr et al. 2010). 

Ultimately the results here rely upon the not unreason- 
able expectation that our knowledge of planetary radii and 
atmospheres will improve dramatically in future years. The 
combination of this improved knowledge and the availabil- 
ity of precision photometry will aid greatly towards breaking 
the degeneracy in these models. 



6 THE UNKNOWN RADIUS AND ALBEDO 

A major source of uncertainty in the discussion thus far 
results from the radius and albedo of the planet. These 
are generally unknown for non-transiting planets and so we 
are using models to describe their dependence on the other 
measured properties, such as mass. Equation [2] shows that 
the planet-to-star flux ratio depends linearly on the albedo 
and quadratically on the planetary radius. Conversely, these 
planetary properties have no impact on either the ellipsoidal 



7 FEASIBILITY DISCUSSION 

The models presented here show that changing the mass of 
the companion will result in unique signatures which can 
be used to constrain the mass and subsequent properties of 
the companion. Such a detection of these signatures present 
a signiflcant challenge to instrumentation requirements and 
our current understanding of the mass-radius relationship. 
The precision requirement for successful detection of the sig- 
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natures for the two examples provided is photometry with 
an accuracy of ~ 10~^. The photometer for the Kepler mis- 
sion is designed to achieve high-precision photometry over 
the 6.5 hour window of a transit, but is not des igned for long- 
term stability over the lifetime of the mission (jBorucki et alj 
|2010| ). The orbits of the radial velocity planets are well un- 
derstood in most cases and so we can accurately predict both 
the amplitude of the predicted phase signature but also the 
phase and times of maximum and minimum flux ratios. In 
contrast, the vast majority of Kepler targets are too faint for 
the acquisition of accurate radial velocity measurments and 
so predicting these times for eccentric orbits amongst Kepler 
targets is more difficult. The afore-mentioned TESS mission 
will provide an opportunity to perform high precision pho- 
tometry on these bright host stars, many of which will have 
known radial velocity planets. The expectation that TESS 
will have sufficient photometric precision to detect down to 
Earth-mass planets around F-G-K stars means that this 
will be very close to the precision requirements for this ex- 
periment. The planned Lyot coronagraph on NIRCam for 
the James Webb Space Telescope (JWST), may be able to 
achieve phase detections for a sample of the most favorable 
targets, though in this case the instrument is optimized to- 
wards young planets around late-type stars. 

Ground-based observations face more challenges in 
terms of correcting for atmospheric stability. However, fu- 
ture generation telescopes will provide opportunities to 
achieve very high precision, such as the European Ex- 
tremely Large Telescope (E-ELT), the Thirty Meter Tele- 
scope (TMT), and the Giant Magellan Telescope (GMT). 
One example of how high precision ca n be achieved from 
such telescopes has been provided by IColon et all (|2010l ') 
whose use of tunable narrow-band filters with the Gran Tele- 
scopio Canarias enabled a photometric precision of < 0.05%. 
Further developments of custom filters, adaptive optics sys- 
tems, and next-generation telescopes will hopefully provide 
a competitive ground-based source for achieving high pho- 
tometric precision. 

A possible alternative approach to detecting the incli- 
nation is astrometric follow-up of th ese targets, such as that 
carrie d out using Hipparcos data bv lReffert fc QuirrenbachI 
l|201ll ). The amplitude of the astrometric signal is roughly 
proportional to the mass o f the secondary component. We 
refer the reader to iKand l|201ll 'l for a more concise com- 
parison of exoplanet detection methods in the long-period 
regime. Thus the examples shown in Section 12.11 will lead 
to a factor of ~ 5 increase in the expected astrometric sig- 
nal. This still results in astrometric signals at the /xarcsec 
level and a challenge to detect from the ground. How- 
ever, the space-based Gaia mission is predicted to have a 
single-measurement as trometric precision of 5-5.5 /.larcsec 
jCasertano et al.]|2008l ) and will therefore contribute greatly 
to siftin g high-mass compan ions from the current exoplanet 
sample (|Sozzetti et al.ll200ll 'l. 



8 CONCLUSIONS 

Perhaps the main ambiguity that is inherent to detections 
of radial velocity exoplanets is the inclination of the orbit 
which produces a measurement of the minimum mass rather 
than the mass itself. The methods described here primarily 



rely on the physical properties of the companion to distin- 
guish between various classes of orbi ts, as oppose d to using 
the orbital properties as proposed bv lBlaclj (| 19971 ). The dif- 
ferent classes of objects which lie along the mass spectrum 
from planets to low-mass stars is becoming better under- 
stood in terms of their atmospheres and photospheric activ- 
ity. However, the change in mass alone is enough to produce 
distinct signatures which can distinguish planets from higher 
mass objects. 

In practice, detection of these signatures is going to be 
difficult to accomplish, even if one can achieve the needed 
precision. Once the mass increases to the point where hy- 
drogen burning can be sustained, the variability effects are 
numerous leading to a veritable plethora of possible phase 
curves for the combined system. Note, for example, that we 
have not considered the changing gravity-darkening of the 
host star due to a variable star-planet distance. The re- 
quired precision is becoming achievable however with cur- 
rent and planned space-based observing platforms. Note 
that, even though the data from the Kepler mission is an 
excellent example of such exquisite precision, the mission is 
searching for transit s of objects around relatively faint stars 
iBorucki et al.ll2010l ). Hence these targets tend to have small 
prospects for a radial velocity orbital solution and no am- 
biguity with regards to their inclination. The kind of post- 
discovery analysis suggested here is therefore more suitable 
towards missions which can easily target the bright stars 
that comprises the bulk of radial velocity exoplanet host 
stars. Missions such as TESS will be able to monitor such 
stars as well as conducting the survey for new transiting 
systems (|Deming et al.ll2009l ). 
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